﻿ ANALIZA SPECTRALĂ A SEMNALELOR ALEATOARE 1 Scopul lucrării Se studiază caracterizarea în domeniul frecvenţă a semnalelor aleatoare de tip „zgomot alb” şi „zgomot roz” şi aplicaţiile acesteia la determinarea modulelor răspunsurilor în frecvenţă ale unor sisteme liniare şi invariante în timp 2 Introducere Modelarea matematică a zgomotelor care apar în dispozitivele şi circuitele electronice, dar şi a semnalelor vehiculate de sistemele de transmisie a informaţiei au necesitat introducerea noţiunii de semnal aleator, echivalentă cu noţiunea de proces aleator sau stochastic din teoria probabilităţilor Pentru a defini un semnal aleator se consideră o experienţă oarecare Prin rezultatul unei experienţe se înţelege una din posibilităţile de realizare a acesteia Mulţimea rezultatelor posibile se va numi în continuare spaţiul eşantioanelor şi va fi notat cu S Un semnal aleator este deci o colecţie de semnale uzuale în timp continuu, numite traiectorii sau realizări Procesele aleatoare sunt semnale, cu două proprietăţi: 1 ele sunt funcţii de timp 2 ele sunt aleatoare, în sensul că înainte de a realiza un experiment, nu este posibil să descriem „exact” forma de undă ce se va genera iul realizărilor Spaţ posibile, S 1s• () 1kxt ) ( 1xt 2s• t s• xt () n ( ) 2xt 2k t () nkxt ) ( nxt t τ 1t kt Figura 1 Un ansamblu de realizări posibile 1 Spaţiul realizărilor posibile conţine, ca puncte realizările procesului, ca funcţii de timp Un astfel de spaţiu sau mulţimea funcţiilor de timp se numeşte proces stochastic sau proces aleator Este evident că se consideră noţiunile de distribuţii în probabilitate ale diferitelor evenimente posibile Evenimentul, sau realizarea, constituie producerea unui anume semnal 3 Procese staţionare iul S i se va asocia o funcţie, cu durata limitată în timp: Fiecărui punct s din spaţ ,Xts, TtT−≤≤ (1) () te şi intervalul de observare Durata 2T se mai numeş Dacă punctul s este fixat, ss= ia de timp (,Xts se mai numeşte şi realizare (de durată ) j, funcţ limitată) sau funcţie eşantion: ,xtXts= (2) () () jj 1, 2, ,xtj n=… Fixând timpul, În figura 1 se arată o mulţime de funcţii eşantion (realizări) ( ) {} j , mulţimea de valori: ktt= …… (3) ()()()()()() {}{} 12 1 2,,, ,,,,,,kk nk k k kkxt xt xt Xts Xts Xts= este o variabilă aleatoare Prin urmare, un proces poate fi privit ca şi o mulţime de variabile ,Xts Pentru simplificarea notaţiilor nu se evidenţiază s şi aleatoare, indexate după timp: ( ) {} procesul se notează simplu cu () Xt Se consideră un proces aleator strict staţionar ( Xt Prin definiţie media procesului ()Xt este ) speranţa matematică a variabilei aleatoare ( Xt şi se notează cu (tμ: ) ) X ∞ (4) == () ()() {} XXttEXt xp xdxμ () ∫ −∞ Xt, pentru t fixat unde p este densitatea de repartiţ ) Xtie a variabilei aleatoare ( () Pentru un proces strict staţionar este valabilă egalitatea: μ= (5) () XXtμ adică media unui astfel de proces este o constantă i fie (,pxx densitatea de Se consideră în continuare două momente fixate, t ş ) 1i 2t ş12,Xt Xt () ( ) 12 i (Xt Atunci media variabilei aleatoare repartiţie comună a variabilelor aleatoare ( )) 1Xt ş2 i (Xt este: produs, ( )) 1Xt ş2 2 ∞∞ (6) ,EXt Xt xxp xx dxdxμ== () ( )() {} () ( )() ( ) 121212 12 1212,,Xt Xt Xt Xt∫∫ −∞ −∞ μ se numeşte funcţie de corelaţie Funcţia care asociază fiecărei perechi ( ) 12,tt valoarea () ( ) 12,XtXt statistică a semnalului aleator şi se notează ( ) 12,XRtt: ,Rttμ= (7) () () ( ) 1212,XXtXt ionar, (,pxx depinde numai de diferenţa Dacă () Xt este un proces aleator strict staţ ) () ( ) 1212,Xt Xt şi nu de valorile absolute ale timpului Prin urmare, avem: 21tt− =−=, t∀, t∀ (8) ()()() 12 12 2 1,XX XRtt Rt t Rτ Proprietăţile funcţiei de autocorelaţie 1 Valoarea medie pătratică a procesului aleator ( este valoarea funcţiei de autocorelaţie Xt ) calculată în origine: 2 (9) 0REXt= ()() {} X 2 Autocorelaţ ia este o funcţie pară: ∈ℜ τ=−, τ (10) ()() XXRRτ 3 Funcţ ia de autocorelaţie ()Rτ are un maxim în origine: X ℜ 0RRτ≤, τ (11) ()() XX∈ Analiza spectrală a semnalelor aleatoare nu se poate face asupra traiectoriilor individuale, pentru că acestea sunt semnale de putere finită, dar această analiză se poate face pe criterii statistice şi energetice Fie, în acest scop, ( ) nxt, cu n fixat, o traiectorie a semnalului aleator Se consideră traiectoria trunchiată: (12) ,2xttT⎧< () n⎪ () ⎨ Txt= 0, in rest⎪ ⎩ transformata sa Fourier Densitatea spectrală medie de putere a acestui semnal, şi ( ) TXω 2 , se obţine împărţind densitatea sa energetică () la durata semnalului: () TXω TSω 21(13) ω= ()() TTSXω T este, pentru ω fixat, o variabilă aleatoare, definită pe câmpul S Notând Se observă că () TSω cu { m⋅ operatorul de mediere şi făcând T→∞ se obţine o funcţie de ω: } 3 21(14) lim limSES EX== ω () ()() {} {} TTωω TT→∞ T →∞ numită densitatea spectrală de putere a semnalului aleator Conform teoremei Wiener-Hincin, funcţ ia de corelaţie statistică şi densitatea spectrală de putere formează o pereche Fourier: F(RS=τω (15) )() {} Fω este o funcţie pară (vezi relaţia 10) şi pozitivă: () 0S≥ω, Rω∈ (16) () , Rω∈ (17) SS=−ωω ()() 4 Procese ergodice Speranţa matematică {E ⋅ este o mediere de ansamblu, pe toate realizările unui proces aleator } ine un alt tip de medie, efectuată „în lungul” procesului, o medie Xt Se mai poate obţ () realizată pe un eşantion al procesului şi efectuată în timp Această mediere poate fi mai uşor implementată, motiv pentru care se doreşte să ştie dacă există vreo legătură între mediile statistice şi mediile temporale Xt, funcţia eşantion (xt, intervalul de observare fiind Se consideră o realizare a procesului ( )) la T Se consideră că (Xt este un proces staţionar Media temporală a lui ()xt este: de la T− ) T 1(18) = ()() xTxtdtμ ∫ 2T T− este o variabilă aleatoare, depinzând de realizarea ()xt curentă şi de Este evident, că ( ) xTμ durata 2T a intervalului de observare Media sa statistică este: TT 11(19) μμ=== ()() {}{} xXXET Extdt dtμ ∫∫ 22TT TT−− adică μ= (20) () {} xXETμ este o estimare neabătută a mediei μ Prin urmare () xTμX sunt satisfăcute două condiţii: Se spune că procesul () Xt este ergodic în medie, dacă ()Tμ tinde spre media statistică, μ, când T→∞, adică: 1 Media temporală xX limTμμ= (21) () xX T→∞ 4 2 Dispersia lui () , considerat ca variabilă aleatoare, tinde spre zero, atunci când xTμ : T→∞ lim var 0Tμ= (22) () {} x T→∞ O altă medie de interes este autocorelaţia Se poate defini o medie temporală pentru estimarea autocorelaţiei: T 1(23) ,RT xt xtdtττ ()()() x=+ ∫ 2T T− ,RTτ este o variabilă aleatoare, dependentă de realizarea (xt şi de lungimea 2T a ()) x intervalului de mediere Se spune, că procesul ( Xt este ergodic în autocorelaţie dacă sunt satisfăcute două condiţii: ) 1 () lim ,RTRττ= () xX T→∞ lim var , 0RTτ= 2 () {} x T→∞ Pentru a avea deci proprietăţile de ergodicitate trebuie ca procesul să fie staţionar Reciproca nu este adevărată 5 Semnale aleatoare în sisteme liniare Se consideră un sistem cu răspunsul la impuls (ht presupus real Dacă la intrarea acestui ) sistem se aduce un semnal aleator staţionar şi ergodic, atunci semnalul de la ieşire este tot un semnal staţionar şi ergodic şi: 2 (24) ωω= () ()() YXSHSω unde () , (Sω sunt densităţile spectrale de putere ale semnalelor de intrare, respectiv de ) XSωY este răspunsul în frecvenţă al sistemului la ieşire, iar () Hω 6 Zgomotul alb Analiza de zgomot a sistemelor de comunicaţii se bazează de obicei, pe o formă de zgomot idealizată, numită zgomot alb, a cărei densitate spectrală de putere este independentă de frecvenţă Adjectivul „alb” se foloseşte în sensul în care se spune că lumina albă conţine în spectrul vizibil componente de diverse culori, cu aceeaşi intensitate O realizare a zgomotului wt şi densitatea spectrală de putere a procesului este: alb se notează cu () oN(25) = () WSω 2 i se exprimă cu: Parametru N se raportează 0 de obicei, etajul de intrare al receptorului ş 5 (26) 0eNkT= unde T fiind temperatura echivalentă e de zgomot a receptorului Funcţia de autocorelaţie a zgomotului alb este: oN(27) δτ= ()() WRτ 2 Se vede că () 0Rτ= pentru 0τ Ca urmare, două eşantioane prelevate din zgomotul alb, W∀≠ indiferent cât de apropiate sunt ele în timp, sunt necorelate Dacă, în plus, zgomotul alb este şi gaussian , eşantioanele sunt şi statistic independente Zgomotul alb ar avea puterea medie infinită şi deci el nu este realizabil fizic Este, mai curând, un concept ce uşurează mult calculele şi conduce la rezultate foarte apropiate de cele din practică Semnalele reale se numesc „colorate” şi au forme diferite pentru funcţiile de corelaţie şi cea de densitate spectrală de putere Cât timp un semnal are o funcţie de corelaţie îngustă, banda de densitate spectrală a puterii este largă şi în acest caz semnalul este mai apropiat de zgomotul alb Dacă funcţia de corelaţie este largă, banda spectrală a semnalului este îngustă şi semnalul este mai apropiat de un semnal periodic (determinist) 7 Desfăşurarea lucrării Să se determine experimental densitatile spectrale de putere ale semnale de la iesirea unui generator de zgomot: zgomot alb si zgomot roz Avem in acest scop: un generator de zgomot, un osciloscop si analizor spectru in timp real (a se vedea anexa pentru descrierea acestuia) Frecventele centrale f ale filtrelor trece banda se citesc de pe panoul frontal al ci analizorului Densitatea spectrala de putere la iesirea fiecarui filtru este 2 VP⎡⎤ i iS= ⎢⎥ HzB i⎣⎦ unde banda de frecvente si factorul de calitate al filtrului sunt: cif B= i Q + 1ci ciff+ Q= − () 12ci ciff+ -12 Obs: Valoarea de referinta este Vref = 1μV , i e Pref = 10 W pentru R=1 Ω: ⎤ dBV⎡ ef⎣ ⎦ 20 [V] 10VV=⋅ ef ref VV WP[ [] ef in] [] dB 10 logP==[] dB 20 logV= ef WP[] [] refVrefV 6 2 VV [] () ef cu R=1 Ω P= R Analiza incepe apasand pe butonul START al analizorului de spectru in timp real; apoi se apasa butonul MOMENT Valorile pentru fiecare canal se citesc in dB; pentru schimbarea 31 canalului se foloseste comutatorul KANAL, de la 2→ Datele se trec intr-un tabel de forma: 22 i cif[Hz] Q iB[Hz] ,ef iV [dB] ,ef iV [V] [V ]iP [V /Hz]iS 2 25 3 31 5 4 40 4 5 31 unde: i = numarul canalului; f= frecventa centrala a filtrului; ci Q = 4,5 (factorul de calitate al FTB); B= fci/Q – largimea de banda; i Pi = puterea; S = densitate spectrala de putere; i Se reprezinta grafic densitatea spectrala de putere a zgomotului alb, respectiv a zgomotului roz 2 – adica [V /Hz]S functie de f[Hz] ici Anexa Analizor de spectru in timp real Schema de principiu a unui analizor de spectru in timp real se prezinta in fig 1 FTB 1 Det Dispozitiv FTB 2 Det de afisare Atenuator de intrare FTB n Det Fig 1 Caracteristicile filtrelor trece-banda FTBi se intersecteaza ca in fig 2 7 1B2BBB i1i+ 22 22 f f f f 1c2cci1ci+ Fig 2 Se cunosc marimile fci, si se stie ca toate filtrele au acelasi factor de calitate: cif (A1) Q= iB Din fig 2 rezulta: + 1iiBB+ (A2) ff−= 1ci ci+ 2 deci + 1ci ciff+ (A3) Q= − () 12ci ciff+ In continuare benzile filtrelor rezulta din (A1) Notand cu xi(t) semnalul de la iesirea filtrului FTBi, si presupunand ca acesta este o tensiune electrica detectoarele ar trebui sa realizeze operatia: /2T 221 (A4) () θ= ,limef i iuxdθ ∫ /2TT− T→∞ In realitate detectoarele realizeaza operatia: /2T 221 θ= (A5) () () ,ef i iut x dθ ∫ /2T− T T fiind o constanta a aparatului Pentru valori suficient de mari ale lui T se obtine 22 ( (A6) ) ,,const ef i ef iutu== Marimea 2 ,ef iU F= i (A7) iB reprezinta o aproximatie a densitatii spectrale de putere a semnalului analizat in banda filtrului FTBi 8